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Abstract. 

We study by means of density-functional calculations the role of lateral surface 
reconstructions in determining the electrical properties of (100) silicon nanowires. 
The different lateral reconstructions are explored by relaxing all the nanowires with 
crystalline bulk silicon structure and all possible ideal facets that correspond to an 
average diameter of 1.5 nm. We show that the reconstruction induces the formation 
of ubiquitous surface states that turn the wires into semi-metallic or metallic. 



1. Introduction 

There has been a growing interest in semiconductor nanowires jllEllS] for their potential 
use in future nanoelectronic apphcations, such as nanocontacts and nanoswitches. 
Sihcon nanowires (SiNWs) are especially attractive for their possible efficient integration 
in conventional Si-based microelectronics. The use of SiNWs as chemical sensors has 
also been demonstrated. SiNW-based sensors for the detection of NH3 jl], of biological 
macromolecules [H] and for the identification of complementary vs. mismatched DNA [0] 
have been reported. Experimentally, SiNWs are grown from a nanocrystal [3 U\ which 
is used as a seed to direct the one- dimensional crystallisation of silicon. The (typically 
gold) nanocluster serves as the critical point for nucleation and enables the addition of 
reactants that allow the growth, determining the nanowire diameter and orientation. 
Ma et al. jHj achieved the thinnest SiNWs reported so far, with diameters as small 
as 1.3 nm for a wire grown along the (110) direction. A few other groups [71 ITU] 
also obtained SiNWs with diameters below 10 nm and growth orientations including 
(100) [7], (110) mini, (111) U and (112) [9J. 

The one- dimensionality is known to induce a gap broadening effect, due to quantum 
confinement, in H-passivated SiNWs [llj. In this paper we explore on theoretical 
grounds the electronic structure of SiNWs when the lateral surface is left free to 
reconstruct. Silicon surfaces are among the most studied systems of the latest decades, 
both theoretically and experimentally, and it is well-known that the reconstruction 
drastically affects their electronic properties. In Section 13.11 we discuss the geometry of 
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Figure 1. Facet arrangement of a (100) SiNW as dictated by the Wulff's rule. The 
formation of {110} facets allows a smoother matching between {100} facets. 

the surface reconstructed SiNWs considered and the imphcation of the Wulff's rule in 
the case nanoscale one- dimensional systems. In Section K12I we analyse the electronic 
structure of the most favoured reconstruction and the localisation of the surface states 
that have been found to form. 

2. Computational methods 

The calculations presented in this paper have been carried out in the framework of 
density-functional theory (DFT). We have used both a numerical atomic orbital J2] 
and a plane- wave [Hj basis set. We have used a double-C polarised basis set jI2I with 
pseudopotentials of the Troullier-Martins [IS] type and a plane-wave energy cutoff of 
20 Ry [Hj with ultrasoft pseudopotentials jTHj. In both cases the exchange-correlation 
energy was calculated according to the Generalised Gradient Approximation ^Hl- The 
wires that we have studied have a diameter of ~ 1.5 nm and a number of atoms 
ranging from 57 to 171, depending on the supercell size and on the adopted shape 
of the unrelaxed wire section. The Brillouin zone has been sampled according to the 
Monkhorst-Pack scheme with a converged grid of 1x1x4, 1x1x6 or 1x1x12 
k-points, depending on the supercell size. 

3. Results and discussion 

3.1. Structural relaxations: checking the Wulff's rule at the nanoscale 

Convincing experimental evidence shows that SiNWs are constructed around a 
crystalline bulk core [HI HE] • The shape of the section of such a structure will depend on 
the way in which the involved vicinal surfaces match and on how abrupt is the transition 
from one to the other. The case of SiNWs grown along the (100) direction is of especial 
interest, because it determines a rather confiictive situation at the edges where two (100) 
surfaces meet with a 90° angle [see Figure [T] and |21^c)]. In such conditions, it is known 
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from Wulff 's theorem JH] that the formation of a facet is to be expected, because it 
would partly release the stress accumulated at the edge. An example of faceting of a 
(100) SiNW |2ni 1^ is shown in Figure ^ As can be seen, the formation of a {110} 
facet allows a much smoother transition between the two vicinal (lOO)-like surfaces. 
However, as discussed previously by Zhang and Yakobson [221, the case of very thin 
one-dimensional structures such as those that we are treating, the predictions of Wulff 's 
rule should be carefully revised. The size of the {110} and of the {100} facets are of 
the same order of magnitude, thus the surface energy associated to the edges can no 
longer be neglected, as assumed in Wulff 's rule. For these reasons, the first part of this 
work has consisted in identifying the favoured faceting arrangement of (100) SiNW of 
nanometric thickness, thus verifying the validity of Wulff 's theorem in these conditions. 

We have considered two different square-section wires [Figure I2ta) and (c)] and 
two Wulff faceting arrangements [Figure |2fb) and (d)]. The square wire of Figure |2ta) 
is made of {110} facets, while the wire in Figure |2fc) has only {100} facets. The 
wires of Figure |2fb) and (d) have been obtained from those of Figure |2f a) and (c), 
respectively, smoothening the corners according to Wulff 's rule prescription. The way 
in which the tension at the edges is released and the way in which the different Si 
surface reconstructions ({100} and {110} facets) compete among them will determine 
the minimum energy structure. Therefore, given the starting configurations of Figure |21 
very different minimum energy structures are expected. 

The results of the relaxations are shown in the section- view of Figure El As can be 
seen, all but one of the wires maintain a high degree of in-plane symmetry. As expected, 
in the SiNW of Figure Efa) {110} facets are dominant. The same relaxation pattern is 
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followed by the wire of FigureEl^b), with the only notable exception of the missing corner 
atoms, whose removal allow a smoother match between the vicinal (llO)-like surfaces. 
In accordance to the prediction of Wulff's rule, the SiNW of Figure Efb) turned out to 
be more stable than that of Figure a) (the difference amounting to ~ 33 meV/atom). 

In the relaxation of the SiNWs of Figure Efc) and (d) {100} facets prevail. It 
is interesting to note that in this case the effects of a Wulff-like profile are more 
evident than in the former case. Here an actual edge-faceting develops, mediating the 
transition between the vicinal (lOO)-like surfaces by {110} facets [see Figure Efd)]. On 
the other hand, when the edges between the {100} facets are not smoothened, the wire 
favours a symmetry breaking, ending up with the butterfiy-sh.a.ped section of FigureEfc). 
Interestingly enough, the butterfly wire of Figure Efc) has turned out to be slightly more 
stable than the wire of Figure Efd), thus confirming that at such small diameters the 
validity of Wulff's rule should be carefully checked. 

Regarding the influence of facet orientation, we have found that wires dominated by 
{100} facets [Figure El^c) and (d)] are more stable than those where {110} facets prevail 
[(Figure ini^a) and (b)]. The cohesive energies are summarised in Tabled For these 
reasons, the candidate structures for SiNWs grown along the {100} direction appear to 
be those in FigureE^c) and (d), favouring the formation of {100} over {110} facets. 

We have found that for the SiNW of Figure Efd) two different reconstructions of 
the {100} facets are possible [21]. Like in the case of the infinite Si(lOO) surface, the 
reconstruction consists in the formation of a sequence of buckling dimers. What differs 
between the two competing geometries that we have obtained is the pattern followed 
by the Si dimers. In one case they determine a trough in the middle of the facet, while 
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Dominant facets Section shape Cohesive energy (eV/atom) 
(110) square 3.886 
(110) Wulff 3.919 
(100) square 4.006 
(100) Wulff 3.989 



Table 1. Cohesive energies per atoms of the SiNW type illustrated in Fig. O The 
SiNWs where {100} facets prevail are approximately 0.1 eV/atom more stable. 




Figure 4. Minimum energy geometries for the facets of (100) butterfly-like SiNW of 
Figure l^ld). 

in the other one every two dimers on one of the two sides is flipped with respect to the 
symmetric behaviour of the trough reconstruction. The relaxed geometry of the facets 
of the butterfly wire is shown in Figure |3 The symmetry breaking is reflected in the 
relaxation and the {100} facets follows two different relaxation patterns. One of them 
- Figure Efa) - is rather symmetric, with all the dimers of one side flipped; the facet of 
Figure IH^b) presents a regular sequence of buckled dimers only in one of the two sides, 
while in the other one every two dimers is missing. 

3.2. Electronic structure 

Nanowires - and SiNWs in particular - are expected to play an important role in 
future molecular electronics applications. Therefore, a thorough understanding of their 
conductive properties is required. In this section we discuss our results of the electronic 
structure of the butterfly wire [FigureEfc)], as well as of the two different reconstructions 
that we have found for the Wulff-like SiNW [Figure Efd)]. 

In Figure and (b) the band structure diagrams corresponding to the two 
competing geometries for the Wulff-like wire are displayed. It is somehow surprising 
to discover that they are rather different: while the trough reconstruction is strongly 
metallic, with four bands crossing the Fermi level, the flipped dimer geometry is only 
semi-metallic, with one band approaching the Fermi energy with a zero derivative at the 
zone boundary. This is quite striking, at first sight, because the differences between the 
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Figure 5. Band structure diagrams of (a) the trough and (b) the dimer flipped 
reconstruction of the SiNW with a Wulff-like section; (c) the butterfly-like SiNW. 



two reconstructions do not seem to be so relevant. However, when considered in more 
detail, one notices the pivotal role that the flipped dimer has in breaking the surface 
Bloch state that would otherwise form. On the contrary, the surface state forms in the 
trough reconstruction where all the Si dimers follow the same buckling pattern. 

The band structure of the butterfly wire is shown in Figure Efc). Also in this 
case the reconstruction results in a metallisation of the wire surface, with one band 
crossing the Fermi level and a semi-metallic band tangent to it at the zone boundary. 
The formation of these metallic and semi-metallic states is directly induced by the 
reconstruction of the surface dimers. Therefore, they are expected to be localised at 
the wire's outer layers. The dimerisation leads to a dangling bond that hybridises with 
the other dimers' dangling bonds along the nanowire in a surface 7r-bond. The surface 
nature of the semi-metallic and metallic states of Figure Etc) is evident in Figure IHl 
where we have plotted the corresponding wave functions. As can be seen, the metallic 
and semi-metallic state are localised, at opposite sides, along one of the (110) diagonal. 
A similar surface localisation, though more symmetric, following the pattern of the 
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Figure 6. Wave function of (a) the semi- metallic and (b) metallic state, respectively 
labelled as (1) and (2) in FigureEtc). 

overall relaxed geometry, is observed for the Wulff-like SiNW of Figure Efd) (not shown 
here; see Reference j2I])- Therefore, conduction in surface- reconstructed SiNWs will be 
almost entirely sustained by the outer layers of the wire, with a negligible penetration 
inside the wire's core. These results, i.e. metallic or semi-metallic nature and surface 
character of the conduction channels, concern all of the (100) wire types that are likely 
to be obtained. As we have discussed in Section EH] (see also Table H]), the difference of 
cohesive energies among the favoured wire geometries [the wire of Figure Efc) and the 
two reconstructions of Figure Efd)] are very small and the selective growth of one or 
another wire seems a difficult task. 

4. Conclusions 

We have shown that in absence of a proper passivation, the lateral surface of SiNWs 
strongly reconstructs, forming series of buckled dimers with a pattern similar to Si (100) 
surfaces. Depending on the reconstruction, surface states that cross the Fermi level can 
form, leading to metallic or semi-metallic SiNWs. Under such circumstances, doping 
is no longer needed to have highly conducting nanowires. The possibility of tailoring 
SiNWs that are conducting without the need of doping is particularly relevant in very 
thin wires where typical doping concentrations require a precision in the fraction of 
impurities per atom extremely difficult to control. 

Available experiments have not thoroughly explored the possibility of surface 
reconstruction or nanowire conductance without doping. However, these results show a 
promising venue of experimental and theoretical research of pure SiNWs. 
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